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ABSTRACT 

We present observational estimates of AF^™ P in a sample of 28 Galactic Globular Clusters (GGCs) 
observed by HST. The photometric accuracy and the sizable number of stars measured in each cluster 
allowed us to single out the RGB Bump both in metal-poor and in metal-rich GGCs. 

Empirical values are compared with homogeneous theoretical predictions which account for both H and 
He burning phases over a wide range of metal abundances (0.0001 < Z < 0.02). We found that, within 
current observational uncertainties on both iron and a-element abundances, theory and observations are 
in very good agreement, provided that the metallicity scale by Carretta & Gratton (1997) as extended by 
Cohen et al. (1999) is adopted. Interestingly enough, we also found that both theoretical and observed 
values show a change in the slope of the AU^g mp -[M/H] relation toward higher metal contents. 

Subject headings: Globular Clusters: general — stars: evolution — stars: horizontal branch — stars: 
Population II 
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1. INTRODUCTION 

The RGB Bump is an intrinsic feature of the RGB Lu- 
minosity Function (LF) of GGCs. It appears as a peak in 
the differential LF or as a change in the slope of the cumu- 
lative LF. The presence of the Bump is due to the fact that 
during the RGB evolution the H-burning shell crosses the 
chemical discontinuity left over by the convective envelope 
soon after the first dredge-up phase. Since its first detec- 
tion in 47 Tuc (King, Da Costa & Demarque 1985), it be- 
came the crossroad of several theoretical and observational 
investigations (Alves & Sarajedini 1999). The detection of 
the RGB Bump was mainly hampered by the size of the 
available samples of RGB stars. This is particularly true 
for the most metal-poor clusters, where the Bump moves 
toward brighter magnitudes and therefore less populated 
RGB regions. Only recently was this feature firmly de- 
tected in a large set of both Galactic (FP; Brocato et al. 
1996) and extragalactic stellar systems (Desidera 1999). 

The analysis by FP showed that the AU^g mp 5 values 
predicted by theory were 0.4 mag brighter than the empir- 
ical estimates in a sample of 11 GGCs. This observable, 
which does not depend on the distance modulus, on the 
reddening, or on the calibration of photometric data, is 
a key parameter for assessing both the accuracy and the 
plausibility of the physical assumptions adopted for con- 
structing the evolutionary models. As a consequence the 
discrepancy found by FP needs to be understood. In order 
to explain this mismatch between theory and observations 
Alongi et al. (1991) suggested the inclusion of convective 
undershooting at the boundary between the convective en- 
velope and the thin H-burning shell, whereas Straniero, 



Chieffi & Salaris (1992) called attention on the role that 
the global metallicity could play on the estimate of this 
observable. The latter hypothesis was confirmed by CS 
who found that standard models agree quite well with the 
observed AV^ mp values in 8 GGCs for which high res- 
olution spectroscopic determinations of both [Fe/H] and 
[a/Fe] were available. The significance of this result was 
limited by the small number of clusters taken into account 
and by the small metallicity range that they cover. 

The main aims of this investigation are to provide new 
homogeneous measurements of both the Bump position 
and the corresponding AV^g mp value in a large sample of 
GGCs covering a wide metallicity range (—2.1 < [Fe/H] 
< —0.2); and also to compare the new measurements 
with the predictions of standard evolutionary models. At 
present, the most homogeneous photometric sample comes 
from HST. Two of us (GP and MZ) are involved in a 
program aimed at collecting F439W and F555W WFPC2 
images for all the GGC with (m — M)b < 18.0 not yet 
observed with HST. Ten of these clusters have been ob- 
served during Cycle 6 (GO6095) and nine more during the 
ongoing Cycle 7 (GP7470). Similar observations for nine 
additional clusters are available on the HST archive. 

In §2 we present this sample of 28 GGCs, and dis- 
cuss the data reduction strategy as well as the approach 
adopted for estimating both the zero age horizontal branch 
(ZAHB) and the RGB Bump visual magnitudes. Section 
3 deals with the comparison between theory and observa- 
tions, while in §4 we briefly summarize the main results. 

2. DATA SELECTION AND REDUCTION 
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zontal Branch (HB) stars located within the RR Lyrae instability strip (see for a detailed discussion Fusi Pecci et al. 1990, hereinafter FP, and 
Cassisi & Salaris 1997, hereinafter CS). 
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The center of all the clusters were observed in the HST B 
(F439W) and V (F555W) bands with the WFPC2. Both 
very short (a few seconds) and long (a few hundreds of 
seconds) exposures are available for both filters allowing 
us to map the evolved regions of the color-magnitude dia- 
gram (CMD). The pre-processing, photometric reduction, 
and calibration to a standard B and V system of the data 
for each cluster were carried out following the same pro- 
cedure as described in details in Piotto et al. (1999a) for 
a subset of images of the GO7470. The photometry for 
all the clusters extends from the tip of the RGB down to 
about 2 magnitudes below the turnoff, and include a num- 
ber of stars ranging from ~ 5000 for the loosest cluster 
to ~ 30000 for the most concentrated ones. A detailed 
description of the CMDs will be presented in Piotto et al. 
(1999b). The large sample of stars, coupled with the high 
photometric accuracy, allows us to identify of the RGB 
Bump even in the most metal-poor clusters. 

Figure 1 shows an example of the Bump identification 
in NGC5824, a metal-poor cluster in our database. The 
upper left panel shows the RGB differential LF. We con- 
structed a histogram with a fixed bin size of 0.15 magni- 
tude (dotted line) as well as the multibin histogram (solid 
line) described in Piotto et al. (1999c). The Bump is 
marked by a vertical arrow. As a further check of the 
accuracy in the Bump identification, the lower left panel 
shows the RGB cumulative LF. The Bump can be easily 
located since it is marked by a change in the slope. The 
right panel shows the CMD of NGC5824 and the arrow 
marks the position of the Bump along the RGB. 

The determination of the ZAHB magnitude (Vzahb) is a 
thorny problem, in particular for those clusters which show 
blue HBs. A standard method to determine the ZAHB 
magnitude for both intermediate and metal-poor clusters 
is to adopt the mean magnitude of the cluster RR Lyrae 
stars. Unfortunately, we could not adopt this method di- 
rectly, since our photometry covers a very short time in- 
terval, and therefore the RR Lyrae were always measured 
at random pulsation phases (Piotto et al. 1999a). 

In order to overcome this problem we have undertaken 
a different approach. For metal-poor and intermediate 
mctallicity clusters we selected from the literature three 
clusters which offer accurate photometry, a large num- 
ber of RR Lyrae stars, and a well defined blue HB tail. 
The template clusters are the following: NGC1851 (Walker 
1998) for clusters in the metallicity range —1.5 <[Fc/H]< 
-1.0, NGC5272 (M3, Buonanno et al. 1994) for clusters 
in the range -1.7 < [Fe/H] < -1.5, and NGC4590 (M68, 
Walker 1994) for more metal-poor clusters. The CMDs of 
the template clusters were artificially shifted in color and 
in magnitude in order to match both the RGB and the blue 
HB tail of each cluster. The mean RR Lyrae magnitude of 
the template cluster was scaled according to the magnitude 
shift adopted to overlap the stellar distributions along the 
HB. The Vzahb magnitudes were estimated by using the 
relation between the mean RR Lyrae magnitude and the 
Vzahb magnitude suggested by CS. It is worth noting that 
this method is totally independent from any zero point dif- 
ference between the template photometry and ours. 

For the most metal-rich clusters ([Fe/H]> —1) with well 



populated red HBs but no RR Lyraes, first we estimated 
the photometric error a v = cn B _ v J2, where o? B _ v ^ is 
the standard deviation of the color distribution of the RGB 
stars at the level of the HB. Then we fixed Vzahb at 3ov 
magnitudes above the lower envelope of HB stellar distri- 
bution. However, in order to provide a consistent deter- 
mination of this parameter over the whole GGCs sample, 
the Vzahb at the level of RR Lyrae instability strip was 
evaluated according to the method suggested by Fullton 
et al. (1995). Figure 1 shows the fit of the same metal- 
poor cluster NGC5824, (full dots) with the template clus- 
ter M68 (crosses). A vertical and a horizontal shift of 2.83 
and 0.08 magnitudes were applied to M68 for matching 
both the RGB and the HB. The M68 RR Lyrae stars are 
plotted as open squares. Taking into account the uncer- 
tainties in the fit, we estimated an error in Vzahb of the 
order of 0.1 mag. 

Table 1 summarizes the cluster observables: column (1) 
gives the NGC and the Messier numbers; column (2) lists 
the visual magnitudes of the RGB Bump and the photo- 
metric error; column (3) gives the mean RR Lyrae visual 
magnitudes estimated according to the method previously 
described; column (4) lists the ZAHB visual magnitude at 
the lower envelope of red HB stars in metal-rich clusters; 
column (5) gives the cluster metallicities according to the 
Carretta & Gratton (1997, hereinafter CG) scale 6 . 



3. COMPARISON BETWEEN THEORY AND OBSERVATIONS 

Figure 2 shows the comparison between theory and ob- 
servations in the [M/H]-AV^g mp plane. The global metal- 
licities were estimated by adopting a mean a enhancement 
equal to 0.3 for clusters with [Fe/H] < -1.0 and to 0.20 for 
more metal-rich clusters. The former value was suggested 
by Carney (1996, hereinafter C96), while the latter, due to 
the paucity of data available in the literature, is a mean be- 
tween the estimates collected by C96 and by Salaris & Cas- 
sisi (1996, hereinafter SC). The observed AV^ mp values 
have been plotted in the top panel according to the CG 
metallicity scale. The AV^ mp error bars have been calcu- 
lated by quadratically combining the errors on Vzahb and 
on Veump- The global metallicity error bars are a lower 
limit of the uncertainties affecting both [Fe/H] and [a/Fe] 
measurements (see C96 and Rutledge, Hesser, & Stetson 
1998, hereinafter RHS). 

The theoretical predictions plotted in the top panel were 
estimated by adopting progenitor masses ranging from 
M/Mq =0.8 to 1.0 and a wide range of global metal- 
licities (-2.3 < [M/H]< 0.0). The initial helium con- 
tents adopted in constructing evolutionary models are the 
following: Y = 0.23 for [M/H]< -0.5,' Y = 0.255 for 
[M/H]= -0.25, and Y = 0.289 for [M/H]= 0.0. Metal- 
poor AV^g mp theoretical estimates up to [M/H]« —0.5 al- 
ready been presented in CS, whereas more metal-rich ones 
were specifically computed, to extend the predictions to 
clusters more metal-rich than 47 Tuc. Basic assumptions 
on the input physics adopted for constructing evolutionary 
models were extensively described in CS and Bono et al. 
(1997 and references therein) and therefore they are not 

6 For those clusters whose metallicity was not provided by CG, the cluster metallicities collected by Harris (1996) were transformed into the 
CG scale by adopting the new calibration provided by Cohen et al. (1999, hereinafter CGBC). In contrast with CG the new scale applies also 
to metal-rich clusters (-2.12 < [Fc/H] < -0.3). 
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discussed here. The reader interested in a detailed dis- 
cussion on the dependence of AV H g mp on mixing-length, 
He content, and element diffusion is referred to CS and 
to Cassisi, Degl'Innocenti, & Salaris (1997, hereinafter 
CDS). Bolomctric magnitudes were transformed into V 
magnitudes by adopting bolomctric corrections provided 
by Castelli, Gratton & Kurucz (1997). 

In order to account for the Veump dependence on clus- 
ter age as suggested by CS and more recently by Alves 
& Sarajedini (1999) in a detailed investigation on HST 
data of 8 SMC clusters, we plotted the theoretical pre- 
dictions for three different ages: 12 (short-dash line), 14 
(solid line), and 16 (long-dash line) Gyr. In the last few 
years a large number of theoretical and observational in- 
vestigations have been devoted to the absolute and the 
relative ages of GGCs as well as to the errors affecting 
such parameters. In fact, they depend on the physi- 
cal assumptions and on the input physics (Cassisi et al. 
1998; Vandenberg, Stetson & Bolte 1996) adopted for 
constructing evolutionary models. As plausible assump- 
tions we adopted an average cluster age of 14 Gyr (Van- 
denberg 1999) and an average uncertainty of ±2 Gyr. 
The change in the slope toward higher metal contents 
shown by theoretical predictions is due to the fact that 
at fixed age an increase in the evolving mass causes a 
smoother increase in the core-mass luminosity relation, 
and in turn in the Veump magnitude. In fact, for metal- 
licitics ranging from Z=0.001 to Z=0.006 the My (Bump) 
and the My(ZAHB) magnitudes changes according to the 
following derivatives: <9My(ZAHB)/d[M/H] w 0.093 and 
<9My(Bump)/<9[M/H] « 0.753, whereas for 0.006 < Z < 
0.02 they change according to w 0.289 and « 1.251 respec- 
tively. 

The data plotted in the top panel of Figure 2 show 
clearly that the discrepancy of 0.4 magnitudes suggested 
by FP is completely removed over the entire metallicity 
range, and indeed the trend of empirical data is well re- 
produced by standard models. This result is even more 
compelling if we take into account that the previous com- 
parisons (FP and CS) were hampered by the small number 
of metal-poor clusters ([Fe/H] < —1.5) for which reliable 
estimates of Veump were available. Figure 2 also extends 
the comparison to [Fe/H]= —0.2; note how the observed 
values in the high metallicity range show the flattening 
predicted by the theory, though with a distribution that 
is somehow flatter. A plausible change in the He content 
(CS) and/or the inclusion of element diffusion (CDS) can 
account for this effect only marginally. As a consequence, 
this result suggests that the metal-rich clusters could be 
younger than the bulk of our clusters (Salaris & Weiss 
1998; Rosenberg et al. 1999). 

In the cluster sample adopted in this investigation there 
are only two clusters -NGC7078 and NGC5694- which are 
marginally in agreement with the theoretical expectations. 
For both clusters the uncertainty on the location of the 
Bump is very small (Table 1). In the case of NGC7078, 
the adopted mean Vrr, magnitude is also in very good 
agreement with the value obtained by Silberman & Smith 
(1995). Therefore for this cluster the discrepance could be 
due to an underestimate of the cluster metallicity and/or 
of the a enhancement (see e.g. CG). An independent esti- 
mate of Vrr for NGC5694 is not available, but any plausi- 



ble assumption on its uncertainty can hardly remove such 
a discrepancy. 

In order to account for the uncertainty on the metal- 
licity scale the middle panel shows the same compari- 
son as the top panel, but the observed points are plot- 
ted according to the Zinn & West (1984, hereinafter ZW) 
metallicity scale. The data plotted in this panel show 
that the AV^ mp values are systematically shifted toward 
lower metal contents when compared with theoretical ob- 
servables. At present, both systematic and observational 
errors affecting the metallicity ranking of GGCs are still 
controversial issues (RHS; C96). This notwithstanding, 
the CG scale is more robust since it relies on recent high 
dispersion spectroscopic measurements and up-to-date at- 
mosphere models. Even though current observational un- 
certainties affects the global metallicity of individual clus- 
ters, data plotted in Figure 2 support the evidence that 
the ZW scale underestimates the cluster metallicity in the 
range —1.7 <[Fe/H]< —1.0. In order to supply a quanti- 
tative estimate of the difference we performed a fit of the 
empirical data (—2.0 <[Fe/H]< —0.5) with predictions at 
14 Gyr. The standard deviation is 0.05 mag for the CGBC 
scale and 0.17 mag for the ZW scale. The latter value is 
almost a factor of two larger than the photometric uncer- 
tainty. 

4. CONCLUSIONS 

We have presented new homogeneous measurements of 
the AV^ mp values for a sample of 28 GGCs observed 
with HST, and a detailed comparison with the theoret- 
ical models. By relying on homogeneous theoretical and 
observational frameworks and on the metallicity scale sug- 
gested by CG and by CGBC we found that, within cur- 
rent uncertainties, observables predicted by standard H- 
and He-burning evolutionary models agree with the em- 
pirical data. This result is further strengthened by the 
fact that this comparison was extended from metal-poor 
to metal-rich clusters (-2.1 < [Fc/H] < -0.2). 

New theoretical predictions for metal-rich clusters show 
a change of the slope of the AV^ mp -[M/H] relation. This 
behavior is supported by the tail of metal-rich clusters in 
our sample and does not depend on the adopted metallic- 
ity scale. Leading physical arguments on the dependence 
of AV^g mp on input physics support the suggestion that 
metal-rich cluster could be younger than the bulk of clus- 
ters in our sample. 

By adopting the ZW metallicity scale we found that em- 
pirical data at low and intermediate metallicity are shifted 
toward lower metallicitics when compared with theory. At 
the same time the agreement between theory and observa- 
tions supports the use of a AV [ f B imp - metallicity relation 
for constraining the cluster metallicity (Desidera 1999). 
However, we note that such a relation relies on the assump- 
tion that all GGCs are coeval within ±1 Gyr (Stetson et 
al. 1999), and that the intrinsic accuracy is of the order of 
0.15 dex provided that the AV^ mp values are measured 
with an accuracy of 0.10 mag. 

We are deeply indebted to E. Carretta and R. Gratton 
for providing us with the extension of the CG metallicity 
scale to [Fe/H]= —0.3 in advance of publication. 
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Table 1 
Cluster Observables 



Object 


Vr 


Vrr 


Vrhb 


\FelH\co 


(1) 


(2) 


(3) 


(4) 


(5) 


104 47Tuc 


14.57±0.02 




14.16 


-0.70° 


362 


15.47±0.02 


15.33 




-1.15° 


1851 


16.16±0.02 


16.05 




-1.14 6 


1904 M79 


16.00±0.04 


16.21 




-1.37° 


2808 


16.31±0.03 


16.25 




-1.24 6 


5634 


17.77±0.03 


17.95 




-1.66 6 


5694 


18.49±0.04 


18.64 




-1.70 6 


5824 


18.10±0.03 


18.46 




-1.69 6 


5927 


17.37±0.04 




16.83 


-0.31 6 


6093 M80 


16.12±0.03 


16.36 




-1.48 6 


6139 


18.30±0.03 


18.40 




-1.50 b 


6205 M13 


14.70±0.04 


14.99 




-1.39 a 


6235 


17.24±0.05 


17.31 




-1.27 6 


6273 M19 


16.77±0.04 


16.91 




-1.53 6 


6284 


17.36±0.05 


17.25 




-1.20 6 


6287 


16.60±0.05 


17.21 




-1.88 6 


6293 


16.04±0.05 


16.47 




-1.76 b 


6342 


18.00±0.05 




17.56 


-0.57 6 


6356 


18.53±0.03 




18.05 


-0.44 6 


6362 


15.60±0.02 


15.35 




-0.96 a 


6388 


17.69±0.04 




17.31 


-0.53 6 


6441 


18.46±0.04 




17.99 


-0.46 6 


6522 


17.01±0.04 


17.06 




-1.38 6 


6624 


16.68±0.02 




16.19 


-0.36 6 


6652 


16.44±0.02 




16.11 


-0.86 6 


6934 


16.85±0.03 


17.01 




-1.40 b 


6981 M72 


17.13±0.04 


17.26 




-1.40 6 


7078 M15 


15.41±0.04 


15.86 




-2.12 a 



a CG metallicities from high resolution spectroscopy. 
b Clustcr metallicities based on CGBC scale. 




Fig. 1. — Example of the empirical measure of the Bump and ZAHB location. The arrows indicate the adopted Bump 
magnitude. Panel a: Differential LF; Panel b: cumulative LF; Panel c: position of the Bump in the CMD; Panel d: match 
between the HBs of NGC5824 and of the template cluster M68. 
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Fig. 2. — Comparison between the theoretical and empirical values of AV^g mp as a function of the global metallicity. Top 
panel: the empirical data plotted according to the CG metallicity scale (see text for more details). Bottom panel: Same 



